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1. Introduction 

Topological defects are interesting classically stable structures which are present in many 
field theories. Their stability has a topological origin. This explains why they are relevant 
for many areas of Physics. One of the most fascinating class of defects is given by Abrikosov 
vortices fl|]. In three-dimensional space they arise as one- dimensional (string-like) struc- 
tures acting as magnetic flux-tubes in type II superconductors. In a simplified setting, they 
are minimum energy configurations in a two-dimensional Abelian-Higgs gauge model 0|. 
Their stability arises from magnetic flux conservation. Apart from Superconductivity in 
Condensed Matter Physics there are many areas of Physics where this type of solutions 
might be relevant. One of the most fascinating is in Cosmology, in which they have been 
studied in connection to Cosmic strings for quite a while |, 1 |. There seems to be 
a wide class of possible extensions of the Standard Model || ^] which possess this kind 
of solutions. Furthermore, other extensions involving extra spatial compact dimensions 
do provide additional situations in which vortex solutions can appear and become rele- 
vant §§, ||, 10 1 . They have also been studied within the standard model itself in connection 



with the dual-superconductivity scenario of Confinement. 

Being interesting mathematical objects in their own right and having shown a factual 
and/or potential interest in Physics, vortex-solutions have been studied for a long time [11]. 
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A recent book [12] compiles many results accumulated over the years. The most fascinating 
situation occurs for a particular value of the Higgs self-coupling marking the transition from 
type-I to type-II superconductors. In this case, known as the Bogomolny limit [13], the 
field equations (expressing the condition of minimal potential energy per unit length of the 
string) reduce to the first order Bogomolny equations. Taubes [14] has demonstrated that 
given q points (the location of the zeros of the Higgs field) there exist a unique solution to 
the Bogomolny equations (modulo gauge transformations). Despite this simplification there 
is no analytical expression for the solutions. For spherically(in 2-d) symmetric solutions, 
including the single vortex case, the asymptotic behaviour and a Taylor expansion in the 
distance to the centre |lj| has been obtained. 

In Ref. ]l(J we proposed a method to construct the solutions on the torus by means of 
an analytic expansion in a parameter measuring the departure of the area from the critical 
value (twice the flux, in units of the Higgs mass). It is known that for smaller areas there 
are no solutions of the Bogomolny equations. The critical value of the area is the only case 
for which there is an explicit analytic solution of the Bogomolny equations, although a fairly 
dull one (vanishing Higgs field and constant magnetic field). In Ref. [16] we computed the 
first 51 terms in the expansion and found the method capable of describing the profile of 
the solution even for fairly large torus areas, and arbitrary positions of the Higgs zeros. 

We should stress here that, as mentioned in Ref. [16], the usefulness of the idea goes 
much beyond the torus case. Indeed, Bradlow [17] discovered that the Bogomolny equa- 
tions can be generalised to a class of gauge-Higgs systems in compact Kahler manifolds. In 
all cases there is a critical volume and the pattern is repeated. This critical case is referred 
as the Bradlow limit, and corresponds to a particular value of a parameter appearing in 
Bradlow's formulation related to the volume. Our approach is that of expanding around 
this value of the Bradlow parameter and is easily generalisable to other Kahler manifolds, 
besides the torus. More general instances of this approach for higher dimensional theories 
are also possible. As a matter of fact, our proposal originated from the study of self-dual 
non-abelian gauge fields on the torus. In that case there are also no known analytic so- 
lutions to the self-duality equations except for constant field strength solutions occurring 
for particular sizes of the torus. In Ref. [18] one of the present authors and collaborators 
showed that one could consistently expand the solution for arbitrary sizes in a power series 
in a parameter measuring the departure from the critical aspect ratios (which allow for 
analytical solutions). The general setting, to be addressed elsewhere fl~9| involves pertur- 
bations of the metric tensor about the particular values allowing for constant (uniform) 
solutions. Finally, we should also comment that our method has also proven useful for the 
study of non-compact manifolds. For example, if one is interested in studying vortices on 
the two-dimensional plane then the torus can be considered as an infrared cut-off to be 
removed in the infinite area limit limit. Indeed, our results of the previous paper and the 
present one show that this methodology is competitive with other alternatives. 

One of the advantages of having an analytic method is that it can be generalised with 
fairly modest effort to other questions involving vortices. The present paper is devoted to 
one of the most interesting applications which is that of studying the scattering of vortices. 
Previously, numerical methods had been employed to study this point [20, 21, 22]. Indeed, 



- 2 - 



what we will do is to make use of our expansion to compute the metric of the manifold 
of solutions of the Bogomolny equations on the torus. This leads, within the geodesic 
approximation 23], to the determination of the classical dynamics of these vortices. The 
next section will be devoted to the presentation of the formalism and derivation of our main 
formulas. In particular, we derive the general structure of the metric tensor to any order, 
as a polynomial in certain coordinates of the moduli space. The explicit computation of 
the first two orders for arbitrary number of vortices are collected in the Appendix. We 
then focused in the particular case of two- vortex scattering. The corresponding results are 
presented in Section ||. In particular, we give the value (up to machine double precision) of 
the coefficients of the polynomial which defines the metric up to fifth order in the expansion. 
This should provide a very good approximation to the metric not too far from the Bradlow 
limit. As we approach the infinite area limit more and more terms of the expansion are 
required. We went up to order 40 in our expansion to compute the salient features of 
the scattering of vortices on the plane. Here, instead of computing the coefficients of the 
polynomials, we defined a grid in the moduli space and computed the metric at these points 
(up to order 40). From these values we determined the metric tensor for vortices in the 
plane with an accuracy better than one per mille. Finally, in the last section of the paper 
we give a summary and conclusions. 



2. Critical vortex dynamics in the geodesic approximation 
2.1 Vortices and the Geodesic approximation 

Let us begin by recalling the Lagrangian density of the abelian Higgs model in four- 
dimensional(4L>) Minkowski space-time: 

C = -\f, v F^ + l -Tp^){D^) - ^(|0| 2 - l) 2 (2.1) 

where is a complex scalar field, and = da — iAn is the covariant derivative with 
respect to A^, the U(l) gauge potential. There exist stationary, z-independent solutions of 
the classical equations of motion characterised by a quantised magnetic flux 2irq (q being 
an integer) through the x-y plane, known as multi-vortex solutions. The Nielsen-Olesen Q 
vortex (N-0 vortex) is the simplest solution, having q = 1. Its potential energy distribution 
has the shape of a thick string parallel to the z-axis having a finite string tension (energy 
per unit length). A particularly interesting case, which we will refer as "critical", occurs 
for A = 1. In that case the equations of motion reduce to first order equations, known as 
Bogomolny equations, that in our units and for positive flux are given by: 

{Di + iD 2 )§ = (2.2a) 

B = \(\-\^) (2.2b) 

where B is the z-component of the magnetic field. The corresponding solutions, for which 
the string tension is now given exactly given by ir\q\, will be referred as critical multi- 
vortex solutions (in the literature they are sometimes called self-dual vortices). The space 
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of multi-vortex solutions is continuous and can be parameterised by 2q real numbers, which 
can be naively interpreted as the the x-y locations of q N-0 critical vortices. The potential 
energy profile for widely separated vortices supports this interpretation. Individual N-0 
vortices are exponentially localised (energy density decreases exponentially away from its 
centre) with a typical size of order 1 in our units. For vortex separations of this order and 
smaller, the field distribution begins to depart sizably from a naive superposition of N-0 
vortices, but still the Higgs field vanishes at q locations in the x-y plane [14]. The location 



of these zeroes will hence serve as coordinates labelling the space of multi-vortices. 

If we are interested in dynamical processes, we can consider solutions of the classical 
equations of motion which tend to well-separated vortices in the past. This can be regarded 
as a scattering process involving two or more vortices. Numerical experiments have been 



performed by integrating the equations of motion to describe this evolution process [22, 24 



21]. For small initial velocities of the vortices it was suggested [23| that the motion can be 
well-approximated by the motion within the space of multi-vortex solutions, as this costs 
no potential energy. This is referred as geodesic motion [E3, EfJ. This is confirmed to be 



rigorous in the limit of small velocities [26], and numerically (21]] it has been shown to be 
a good approximation even up to relatively high velocities (0.4c). 

Within the context of the geodesic approximation, the kinetic term in the action be- 
comes a quadratic form in the time derivatives of the position of the q zeros of the Higgs 
field. From a geometrical stand-point this quadratic form can be regarded as a metric in 
multi-vortex space. The trajectories which are solutions of the equations of motion are the 
geodesies corresponding to this metric. This explains the name given to this approximation. 

Notice that the non-trivial scattering of vortices destroys the naive interpretation of 
vortices as non-interacting objects. The latter image is based on the fact that the potential 
energy is independent of the position of the vortices. Instead, a very nice physical interpre- 
tation of vortices can be given p7fl in which a vortex is considered a point particle with a 
magnetic charge and a magnetic dipole moment. With this in mind, the interaction of two 
static vortices can be shown to be zero, but when the vortices move along some trajectory, 
it can be shown that the interaction between them is exactly what is expected for large 
separation between vortex centres. 

Despite the nice formalism and results described in the previous paragraphs, the main 
difficulty in studying both static and dynamical properties of vortices is that there are 
no known analytical solutions of the Bogomolny equations, not even for the q = 1 case. 
Most studies rely on numerical methods which often make use of symmetries to reduce 
complexity [^, [0| [2|] . 

2.2 Vortices in compact manifolds 

One can extend the previous results to more general space-times. For example, one can 



generalise the Bogomolny equations [17] to a class of systems with Higgs and gauge fields 
in spaces which are arbitrary Kahler manifolds. In this construction there is a natural 
parameter, which we will call the Bradlow parameter, which measures the spatial volume 
in units of the inverse photon or Higgs mass (which are equal at criticality) . In what follows, 
we will restrict to 2 dimensional compact spaces, extending the previously described case 
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of Euclidean space. These solutions can be regarded as solutions in higher dimensional 
spaces, not depending on the additional coordinates. The space of solutions (moduli space) 
of the corresponding Bogomolny equations turn out to be 2g-dimensional manifolds. The 
kinetic term induces a metric and, hence, a Riemannian structure on them. One can 
also introduce a complex structure and show pi] ] that they become Kahler manifolds. 
Some global properties, such as the total volume of these manifolds of moduli p9| , can be 
determined in terms of the properties of the 2-d space manifold. 

Obviously, the simplest (and usually the most interesting) cases arise for the sphere 
S 2 and the 2-torus T 2 , which have been explicitly addressed in the literature |3(], gl], p2|] . 
The interest of these extensions goes beyond the simply academic or purely mathematical 
one. The torus might, for example, represent the behaviour of vortices on a dense envi- 
ronment. Furthermore, there is an increasing interest in vortex solutions within physical 
scenarios with extra dimensions || ||, [llj . Frequently the latter are compact and have two- 
dimensional cycles, so that these types of solutions are present. A final motivation could 
follow from our work, which shows that one can use the torus results to obtain information 
about vortices on the plane. From now on we will concentrate on the flat torus case. 

Let us consider a two dimensional torus, described as a single patch with coordinates 
x\ and X2 having orthogonal periods of lengths l\ and I2 respectively. We consider the 
metric to be Euclidean (other constant metrics can be easily considered), so that the total 
area is given by A = l\l 2 - We can also introduce the complex coordinate z and its complex 
conjugate z (in general, we will denote complex conjugation with an "overline" symbol), 
and their corresponding partial derivatives: 

XI + IX 2 _ Xx- IX 2 , n Q v 

z = - ; z = (2.3a) 

d = d 1 -id 2 ; d = d 1 +id 2 (2.3b) 

Now we will briefly revise our formulation, as given in Ref. Jlq] . The vector potential (as 
any one-form in this complex manifold) can be combined into a complex valued function 
((l,0)-form) A(z,z) = A\(xi) — zA2(xi). However, this is not necessarily periodic, as 
a one-form should. Its periodicity properties depend on the transition functions on the 



bundle (See Ref. [33] for a simple introduction to the subject of gauge fields on the torus). 
However, if A^ ' is a specific gauge field satisfying the boundary conditions, then A — A^ 
is periodic and we can make use of the Hodge decomposition theorem to parametrise our 
gauge field as: 

A = -idH + v + A^ (2.4) 

where H = h — ig is a complex- valued periodic function of the variables z, z, and v is 
a complex constant labelling the space of harmonic forms on the torus. The function 
g describes periodic gauge transformations which are continuously connected with the 
identity. The remaining discrete set of gauge transformations imply that one can regard v 
as the complex coordinate of the dual-torus (having lengths j^, -^). The periodicity in v 
can be easily guessed by computing the corresponding Polyakov loops (the holonomy). 
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Although a good deal of our formulas do not require a choice of the transition functions, 
for actual calculations one is required. In that case we will take the following quasi- 
periodicity conditions for the field (section of the bundle): 



4>(xi + h,x 2 ) 



i-nq- 



-i-nq-r 

e ! i 



2 <p[Xx,X<2) 

i 



(2.5) 



where q is the first Chern-number of the bundle. 

There is a natural choice for which we will adopt, namely that which leads to a 
uniform field strength, = f = 2*2. With our choice of the transition functions A^ ' is 
given by 

= -tfz (2.6) 



We will complement our parametrisation of the gauge field Eq. 2A , with the following one 
for the Higgs field: 

= V~ee~ H X (2.7) 

where e = 1 — 2/. The function x satisfies the same boundary conditions as the Higgs field 
(Eq. |2.5| for our choice of transition functions). 

Finally, we can write down the Bogomolny equations in our parametrisation: 



(d + fz)x = ivx 



(2.8 



ddh=-[\ 



e~ 2h \x\ 2 



(2.8b) 



These equations can be solved sequentially. First, one solves for x using the first equation. 
Plugging the result into the second equation (known as the vortex equation) one solves for 
h. As we will show, the solution of the second equation is unique given x- Thus, the whole 
structure of the moduli space of multi-vortex solutions resides in the first equation. 

The constant v can be taken as one of the coordinates of the moduli space. The 
remaining coordinates can be associated to the space of solutions at v = 0. This is so 
because given a solution x of Eq. 2.8a for some value of v ^ one can associate to it a 
unique solution x^ f° r v = as follows 



X {0) (z,z)=T v X (v) (z,z) 



IV _ IV 



(2. 



Which is a translation up to a gauge transformation and a rescaling by a constant factor. 
This result extends to Eq. |2.8b by also translating the function h, and adding a constant 
to it. This shows that the moduli space of solutions has the structure of a fibre bundle 
whose base space is a torus parametrised by v. 



As explained in Ref. [16|, for v = Eq. 2.8a has \q\ linearly independent solutions, 
which we will label Xi°\ x )- I n that paper several choices of Xi°\ x ) are given. In particular, 
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we will employ the following one 1 : 



Y (0) , 

As+1 



2qk 
h 



1/4 



s/q 




2qirz I2 



where s = 0, 



, q — 1. They satisfy the following orthonormality conditions 

Si 



"7 / & x xf\x)xf ) (?) 



'K> 



(2.11) 



(2.12) 



Using the operator (T v ) 1 one obtains an associated basis of the space of solutions Eq. 2.8a 



for any v within one patch of the dual torus. A general solution of Eq. 2.8a is then given 
by 



X 



» 



(2.13) 



2 CiX i 
i=i 

where q are |g| complex constants. In principle, one might be led to conclude that the q 
and its complex conjugates, together with the coordinates (v,v) provide coordinates of the 
moduli space of solutions. However, notice that the Higgs and gauge field are invariant (up 
to a global gauge transformation) under the replacement: 



KC; : 



h — ► h + log(\n\) 



(2.14) 



where n is an arbitrary complex number. Hence, for fixed v the space of solutions becomes 
topologically equivalent to the complex projective space CP 9-1 and the coefficients q can 
be regarded as homogeneous coordinates in this manifold. This completes the structure 
of the moduli space as a fibre bundle with fibre T = CP y_1 . This result is a particular 
case of the results of Ref. [31] where this structure is proved for vortices in more general 
manifolds. This parametrisation gives us the first set of coordinates that we will use 



Bi = {v,v,Ci,Ci} (i = l,...,q) 



(2.15) 



A standard choice to obtain inhomogeneous coordinates is to divide all the q by one of 
them (say c\). This gives good coordinates in the patches where c\ 7^ 0. For the q = 2 



case, since 



S , we will also make use of the standard spherical coordinates. 



Another quite natural choice of coordinates of the moduli space is given by the position 
of the zeros of the Higgs field {z r }. In Ref. [16], we verified that the function must 



have \q\ zeros (counted with multiplicity), and that the centre of mass of these zeros is 
located in the point zjP^ = j(h + ifa)- Since the application of the operator T v does not 
change the number of zeros, the function w iH also possess \q\ zeros obeying 



E 



V 



(2.16) 



lr rhe symbol ■& 



(z\t) represent the theta function with characteristics (a, b), and is given by 



\ \ ^ Z7TT 



(2.10) 
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Hence, the coordinate v of the moduli space can be seen to be directly given by the location 
of the centre of mass of the zeros of the Higgs field (note that the coordinate v must lie 
in its proper range. This can always be achieved using the fact that the position of the 
zeros are defined up to translations by one period). Then, one can take as coordinates of 
the fibre the position of the zeros relative to the centre of mass, cjj = Z{ — Z c (where 
Z c = ~ z%j- With an appropriate ordering of these relative coordinates, we can use the 
set 

B 2 = {v,v,Ui,uJi} (i = l,q-l) ( 2 -17) 

as coordinates of the moduli space. 

From this second construction, it is clear that this moduli space is a complex manifold 
having v and uj{ as possible complex coordinates. Instead, one might use v and Cj. The 
corresponding change of coordinates (see Ref. for details) is holomorphic. 



2.3 Bradlow parameter expansion 

The main difficulty in obtaining the explicit form of the solutions is the non-linear character 



of the vortex equation Eq. |2.8q . In our previous paper we gave a recipe for solving it by 
an expansion in the parameter e, directly related to the Bradlow parameter. For e = 
the solution is given by a constant value of h. This constant is fixed uniquely by the flux 
condition and the normalisation of \- I n particular, it vanishes if x is chosen of unit norm. 
In this limit the Higgs field vanishes and the magnetic field is uniform. For small positive 
values of e, one expects a solution which departs only slightly from a constant. In general, 
one can write the solution of the vortex equation as a power series in e: 



h = J2h^e n (2.18) 



n=l 



Plugging this into Eq. 2.8b and matching equal powers of e on both sides, one obtains a 
series of linear non-homogeneous equations which allow the determination of the functions 
iteratively. In practice, for the torus case the solution of the linear equations can be 
accomplished by expanding h^ n \x) in Fourier series, whose coefficients can be determined 
by the iterative procedure mentioned previously. The main difficulty is that to obtain these 
coefficients one has to perform a increasingly large number of convolutions of the previously 
determined M m ) (x) . Fortunately, the Fourier modes of these functions fall very rapidly and 
it is possible to compute the dominant ones with 10-15 significant digits by performing these 
convolutions numerically (truncating the infinite sums down to finite sums). 



In Ref. [16] we were able to carry the e-expansion up to order 51 with a fairly modest 
computational effort. Having those many terms in the expansion allowed us to explore the 
convergence of the series. It is, of course, impossible in this way to prove convergence of 
the series. However, even if the series is asymptotic, it is possible to test what precision 
that can be attained for a particular value of e. Here it is important to notice that by 
virtue of its definition, the range of e values that we are interested in goes from to 
1. The upper edge of the interval e = 1, corresponds to the infinite volume limit of the 
torus: the plane. The results of our analysis turned out to be very positive. Meaningful 
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results, comparing well in value and precision with that of other authors, were obtained 
for the worst possible case of the plane (e = 1). Our method, in contrast to these other 
results, gives us the solutions for all torus sizes and can be applied irrespective of the 
value of q or the symmetries invoked. Having an analytical expansion for the solutions has 
many advantages over purely numerical methods. In particular, many different theoretical 
questions involving vortices can be attacked in the same way. This paper develops one 
of these applications, the determination of the dynamics of vortices within the geodesic 
approximation and our Bradlow parameter expansion. In concrete, the goal is to compute 
the metric in the moduli space of multi-vortices on the torus. This will be explained in the 
next subsection. 



2.4 Derivation of the metric 

We will now allow our fields Ai and <j> to depend on time (while remaining z-independent), 
but we will impose that at any time they correspond to a configuration of minimal energy 
within the corresponding flux sector(2-7rg), i.e. a solution of the Bogomolny equations. 
Therefore, the motion is indeed within the moduli space of solutions of the Bogomolny 
equations. 

To study dynamics we have to compute the kinetic energy, which in the Aq = gauge, 
takes the form: 

T= l -j <Px{A 2 + \^) (2.19) 

The potential energy is fixed to be tt times the total magnetic flux, and is time- independent. 
Thus, the dynamics of the system is given by this kinetic term. This amounts to a finite di- 
mensional classical system, which corresponds to geodesic motion within the 2q-dimensional 
manifold of solutions of the Bogomolny equations. To determine the metric, we simply have 
to parametrise Ai and 4> an d their time derivatives in terms of whatever coordinates of this 
manifold we choose, and plug them into the expression for T. 

We will be using the parametrisation of the Higgs and gauge fields introduced pre- 
viously. In addition to the Bogomolny equations, we also need the dynamical equations 
involving time derivatives. Some can be deduced by differentiating the Bogomolny equa- 



tions |2.8| with respect to time. We have 

(d + fz)x {v) =i(vx {v) +vx {v) ) (2.20a) 

ddh = h\<t>\ 2 - e -e- 2h ^ t \ X ^\ 2 (2.20b) 

One has then to add the Gauss constraint (one of Maxwell equations for the system) 

diAi = Im(^) (2.21) 

This equation can be combined with the time derivative of the second Bogomolny equation 
to give a complex equation: 

d~X = -i~4>4> (2.22) 
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In summary, the only equations that remain to be solved are: 



ddh=-[l 
2 



(2.23a) 



OH = ee- Ih x (v) X (v) (2.23b) 
dd + \(f)\ 2 . These equations determine both h and H in terms of x an d 



Where e 

X • Integrating both sides of Eqs. 2,23 , we obtain 



I 



d 2 xe- 2h X ^X {v) 



d z xe- zn x (v> H X 



(2.24a) 



(2.24b) 



that are basic equations for the constant pieces of h and H. Eqs. 2.22 and 2.24 uniquely 
determine h and H given x an d X ■ Hence, parametrising x m terms of certain 
coordinates of the moduli space, one can then use Eqs. 2.23 , 2.24 , to express h and H in 
terms of these coordinates and their time derivatives. 

Using our parametrisation, we might express the kinetic term as 



2 



d 2 x 



\v\ 2 + ee- 2h 



v\ + ee" 



HOH 



(2.25) 



This is the crucial formula from which the metric and its perturbative expansion can be 
obtained. It is not hard to show that T is invariant under the replacement of x^ by 
j/( v ) — \x^ v ' for any constant A. This expresses the irrelevance of the normalisation of X ■ 
In Eq. p. 25 the time derivatives of the coordinates appear in v and x^ ■ These two terms 
are not independent variations within our fibre bundle choice of coordinates. One should 
rather use v and x^ ■ For that purpose we need the formula 



d,t 



(%r 1 ix (0) -ff(d-fz+ l v)x {0) ] 



(2.26) 



Using this expression we might also relate the solution of Eq. 2.23b| for v ^ 0, which we 
will note H v , with the corresponding one for v = 0: 



iv _ 
Z+ 2f>* 



iv 

V 



H(z,z) 



IV 

2/ 



(8(z, z) + iv) 



where 5 is the solution of the equation 



OS 



ee 



" 2h x {0) (d-fz) x 



(0) 



(2.27) 



(2.28) 



Now we should plug these relations into the formula for the kinetic term T, and after 
changing variables in the integrals and decomposing x^ m the aforementioned basis, we 
arrive at an expression of the form 



7T 



'g vv (ci)vv + gi V (ci)ciV + g V i(ci)vci + gij(ci)diCj) 



(2.29) 
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Where g a p is the induced metric in the moduli space (we follow the normalisation of 
Ref. |3^| ). Furthermore, gi v and g„ vanish. To show this we focus on the term proportional 
to v in the expression of T. It has the form 



(2.30) 



where 



(<9 - fz)x 



(0) 



x (0) 



X 



(0)h 



ce 



-2h 



-dddH + edH - 2d(dH dh) 



(2.31) 



This, being the derivative of a periodic function, has a vanishing integral. By hermiticity 
the same holds for the term proportional to v. In a similar fashion (by partial integration) , 
one can show that the term proportional to \v\ 2 is a constant, to be given below. Finally, 
we arrive to the following expression for the metric 

A 2 



ds z = " * dvdv H — 
A-K z q it 



d 2 xe~ 2h 



\X 



(0)1 2 



O" 



e- 2 ¥ 0) X (0) )} (2-32) 



The first term is the metric associated to the dual torus, while the second is the metric of 
the fibre T . This factorised structure was previously found in Ref. |3l| . 

From now on we might focus upon the metric of the fibre which can be written as 



ds 2 = - 

7T 



d 2 xee- 2h x i0) (x {0) 



I'M = zl 

7T 



d 2 xd 



'x (0) --\ 



X 



(0) 



(2.33) 



Notice that, this time, the integral of a partial derivative does not vanish, because the 
quantity inside parenthesis has singularities (poles) at the zeroes of x . Hence, the metric 
is expressed by the value of this function in the vicinity of these zeroes. This property was 



discovered in Ref. [21] for vortices on the plane. In particular, the metric is twice the sum 



of the residues of the quantity inside brackets of the last equality at its poles. 

Our method of computation of the metric follows from substituting into Eq. 2.33 the 
expression of h and H as a power expansion in powers of e. This expansion follows from 



solving Eqs. |2.23; and |2.24| order by order in e. The resulting expression of the metric in 
the fibre T has the form 



9ij(c k ) = Yl 9ij(c,c)t N 



(2.34) 



N=0 



The computation of the first two orders is detailed in Appendix [A]. The first non-zero term 
in this series, is given by 



91 



7T C 



(2.35) 



where Nell 2 



This is just a multiple of the Fubini-Study metric. This result 



agrees with the metric obtained in this limit in Ref. |3C], but for vortices living in the two 
sphere S 2 . Using this result we can compute the volume of the fibre. 



Vol(^) 



9-1 



M) 
(q-l)\ 



(2.36) 
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This value times the volume of the base-space dual torus matches the calculation done in 



Ref. 1 32], and generalised in Ref. [29] to other manifolds. Furthermore, the result is valid 
to all orders in e. 

In Appendix [A| the second order is explicitly computed. The result involves certain 
constants which can be expressed as an infinite sum over two integers, or as integrals of 
theta functions with characteristics. The general structure of the metric to all orders will 
be described in the next subsection. It includes certain constants given by an increasing 
number of infinite sums over integers. Fortunately, these sums are very fastly convergent 
and a few terms suffice to obtain a double-precision (~ 16 decimal places) numerical de- 
termination. This fact enabled our study of the q = 2 case up to fortieth order, to be 



presented in Section 



2.5 Symmetries and general properties of the metric. 

Here we will focus upon the metric on the fibre T expressed in homogeneous coordinates 
Cj. Some of its properties, which were known before, hold order by order in our expansion. 
For example, we know that the metric is Kahler |H]]. As mentioned previously, we also 
know the total volume of the fibre ]p9f| . Since the volume is already obtained at leading 
order, this tells us that the higher order contributions do not contribute to it. 

Another property of the functions gij(c, c) has to do with the homogeneous nature 
of the coordinates, and has been mentioned previously. Translating the property to the 
expansion coefficients we have: 

5>$°=° (2.37) 

3 

AN) 



A more specific property of the terms g\j follows from the e-expansion of Eqs. 2.23 and 
2.24 : Each new order in e is accompanied in the right-hand side by a power of CiCj/\ |c| | 2 . 



Analysing in detail the different equations, we conclude that Hc]] 2 ^^ is a polynomial of 
degree 2N in the variables and its conjugates. This property is quite crucial for our 
expansion, since it allows the exact determination of the metric to any finite order in terms 
of a finite (and small for N small) number of coefficients. Combining this fact with the 
previous properties, a reduction in the number of free coefficients follows. For example, 



to leading order we have a polynomial of degree 2 in A. Property Eq. 2.37, reduces 
the indetermination to a single overall pre-factor, which is fixed by the volume. Similar 
simplifications occur for the second order result, described in the Appendix. 

A further reduction in the number of free coefficients follows from symmetries. These 
are isometries, which have an origin in a geometrical symmetry of the problem: transfor- 
mations that leave the torus invariant. This includes reflections and rotations by ir. If the 
shape of the torus is square 1% = I2, we also have rotations by \. In addition, we have 
translations. This might look surprising, since we are restricting to a fibre which has a 
fixed position of the centre of mass of the vortices. However, translations by a fraction 
1 jq of the total length in either direction do not change the position of the centre of mass. 
This is a finite group with q 2 elements. 
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In order to exploit the consequences of the symmetries one has to make a choice of the 
basis ■ Here, we will give the results in the basis X s +i = ^Os (defined in Eq. 2,11| ). The 



explicit properties of these functions needed to express the requirements of the symmetries 
are explained in the Appendix of Ref. [16|. Some subtlety is necessary. This has to 



do with the fact that transition functions are, in some cases, not invariant under the 
Euclidean transformations that are candidate symmetries. To recover the symmetry one 
has to accompany the geometrical transformation with a gauge transformation, designed 
to leave the transition functions invariant. We will now present the results. The proofs are 
simple (after consulting Ref. fill ), and are left to the reader. 

The invariance under translations along the x axis, implies that the metric is invariant 
under the replacement 

c s — ► e- 2nis / q Cs (2.38) 
Translations along the y axis are associated with the transformation 

c s — ► Cs+i (2.39) 

where c q = cq. Reflections with respect to the x\ = line are associated with the trans- 
formation 

c s — > c s (2.40) 
Finally the it/ 2 rotations are associated with the transformations: 

9-1 1 

c s — ^c 2 ™ s 'l*cs> (2.41) 
,'=0 v 1 

The square of this transformation provides the rotation by 180 degrees. In the next section 
we will make use of these symmetries for the q = 2 case. 



3. Two-vortex dynamics 

Here we will study with higher detail the specially interesting case of two vortex dynamics. 
As mentioned previously, the fibre T of the moduli space is diffeomorphic to the two-sphere. 
In the first part of this section we will explain the different coordinates that we will be 
using. Then we will present our results focusing on two different regimes. For volumes not 
too far from the Bradlow limit, a few orders in the expansion provide a good description of 
the metric and dynamics. It turns out that the metric is fully determined by fixing a small 
number of real constants. After explaining this general structure, we give the numerical 
value of the 17 constants necessary to fix the metric up to fifth order. With this information 
the corresponding two-vortex dynamics is analysed. On the opposite extreme we have the 
case of very large torus sizes. In this regime, we expect our results to match with the 
case of vortices in the plane. By analysing the results up to 40th order in the expansion 
for a 40 x 40 grid of points in the moduli space we verify that this is the case. Our 
method proves to be strongly competitive in extracting information about the metric and 
scattering of vortices on the plane. 
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Our formalism, given in the previous section, expressed the metric in homogeneous 
coordinates: 

ds 2 = gijdciidcj (3-1) 

and deduced certain properties for the hermitian matrix gij, and its expansion in powers 
of e. In particular, the metric function satisfies c\gn + C2g%2 = 0. Going over to the 
non- homogeneous complex coordinate Z = we get 



„ 2 _ m z) - 

(1 + |Z| 2 ) 2 



dg2 = J^,~J - 



where f(Z,Z) = \\c\\ 2 Ti(g) = ||c|| 2 (gn +522)- It is also useful to express the metric in 
spherical coordinates (6, ip), defined by the following change of variables: 

sin 6* cos ip (3.3a) 



1 + \Z\ 2 
-2b 



cos 9 (3.3b) 



1 + |Z| 2 

where a, and b are the real and imaginary parts of Z [Z = a + ib). In these coordinates 
the metric takes the form: 

ds 2 = {dO 2 + sin 2 6 d^ 2 ) (3.4) 

where, the function f s is the same as / but expressed in the new variables. 

We will also use the more conventional coordinates of the moduli space: the zeros of 
the Higgs field. Let z = ube the location of one of the zeroes (the position of the other zero 
is given by the centre of mass condition). In order for the change of variables to remain 
e-independent, we will consider first w, the position of the zero relative to the torus, given 
by: 

4"7T . . 

w = —u (3.5) 

Now, we can make use of the expressions given in Ref. || to give the holomorphic change 
of coordinates: 

3 (w\i2t) ^qq(w) 
Z(w) = — ; r = — 7 — r 3.0) 

where "di{z\t) are the i th classical Jacobi theta functions 2 and r = £2/^1- F° r ease °f notation 
in the next formulas we will omit the second argument of theta functions, assuming it is 



2 A11 of the four classical Jacobi theta functions can be defined in terms of the theta function with 
characteristics (see footnote [l], page |^) 



3 (*|t) =■& 
®i{z\t) = -$ 






1/2 
1/2 



(z\t) (3.7a) 
(z\t) (3.7b) 
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equal to i2t). Finally, we arrive to the form of the metric in the relative position coordinate 
w: 

(|i? 2 (w)r + |i?3(w)r' 



where / w is obtained from / by change of variables. From here, one can trivially re-express 
the metric using u and u as coordinates. 

After the presentation of the coordinates we will now deduce the implications for the 
function / of the properties of the metric studied in the previous section (The Kahler 
property is obvious in our case). Let us examine how the symmetry transformations act 
on the variable Z. Translations along both axis lead to the transformations Z — ► — Z and 
Z — > 1/Z. Reflection invariance leads to the transformation property Z — ► Z. Hence, 
the function f(Z, Z) must be even and symmetric under the exchange of its arguments and 
satisfy 

f(Z,Z) = f(±, |) (3.9) 
Finally, the transformation associated with tt/2 rotations (if l\ = I2) corresponds to 

" ' " 1 (3.10) 



l + Z 

implying that / must be invariant under this change of variables. 
3.1 Near the Bradlow limit 

Here we will examine the first few orders of the expansion of / in powers of e. These 
terms will describe quite accurately the structure of the metric for small values of e or, 
equivalently, for torus sizes which are not too large compared with the critical area of 
Bradlow. 

As explained in the previous section, the contribution of order to the metric tensor 
I |c| \ 2 g\f^ is a polynomial of degree 2N in q/| |c| |. Going over to the variable Z, we conclude 
that the Nth order contribution to f(Z, Z) is given by 

f (N) (ZZ)- P n{Z,Z) 

J lZ ' Zj " (l + | z |2)(^i) V- LL ) 

where P^(Z,Z) is a polynomial of degree 2N — 2 in its arguments. The symmetries 
imply certain restrictions on the coefficients of this polynomial. Writing Z = a + ib, we 
conclude that P^(Z, Z) is a real polynomial in a 2 and b 2 . This brings down the number of 
coefficients to N{N + l)/2. Further restrictions follow from the symmetry Z — * 1/Z. It is 



easier to express them in terms of the variables 9 and ip defined in Eq. |3.3| . One concludes 



■& 4 (z\t) = I? 



1/2 





1/2 

See Ref. [fi4l for more detailed information about Jacobi theta functions 



(z\t) (3.7c) 
(z\t) (3.7d) 
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0(3) 




0(2) 


4(3) 
^00 


1.57 24 57 43 63 76 56(1) x KT 1 


A (2) 

^00 


2.09 85 65 63 87 72 34 78 07 x 10" 1 


A (3) 


-7.35 97 58 95 32 71 25 x KT 1 


A {2) 


-6.29 56 96 91 63 17 0439 78 x 10" 1 


A3) 
^20 


4.40 39 77 74 02 36 00 x 10" 1 




0(5) 




0(4) 


^00 


3.24 86 80 02 33 12 76 x 10" 1 


^00 


2.56 44 31 24 38 37 31 64(1) x 10" 1 




-1.59 28 19 34 17 36 19 6 




-1.09 63 43 50 68 01 26(2) 


^20 


2.35 30 76 79 79 16 14 




1.12 62 17 59 37 90 85(1) 


^30 


-1.33 17 87 41 82 85 50 3 




-3.23 47 12 72 58 54 03(4) x 10" 1 




2.44 37 72 51 19 98 28 x 10" 1 


^Ol 


-2.62 60 73 46 53 83 54(1) x lO" 1 


^01 


-4.46 53 65 89 63 78 00 x 10" 1 






A (5) 


4.62 92 29 53 30 50 51 x 10" 1 



Table 1: Numerical coefficients which determine the two- vortex metric up to fifth order for l\ = 1%. 



that fs N \o,ip) is a polynomial in cos 2 (#) and sin 2 (9) sin 2 (ip) of degree (N — 1). This is 
exemplified by the first two orders, explicitly calculated in Appendix [A|. The leading order 
term reduces in our case to the a constant curvature metric on the sphere of radius \J~^- 
The next order is given by a sum of three pieces, which are related by the condition that 
the total contribution to the area is zero. 

We will now restrict to the most symmetric case l\ = li- The additional symmetry 
resulting from 90 degree rotations, implies the following general form: 

. j+2fc<iV-l 

/(AO [QM = - y 4f cos 2 ^ 6 sin 4fc cos 2fc 2y (3.12) 

j,k=0 

where A^/ are real coefficients. It is relatively easy to compute these coefficients numer- 
ically up to machine double precision. The results up to fifth order of the expansion are 
collected in table |l[ 

From Eq. 3.12j one sees that the number of terms grows like iV 2 /4, although from the 



table we see that some coefficients turn out to be 0. In addition, there are relations among 

(2) (2) 

them. For example, A\q = —3A 00 . This is a particular case of the restriction following 
from the vanishing of the volume contribution beyond the leading order. This leads to the 
following relations among the coefficients: 

V A<"> T{J + ^ [m]2 - (3 13) 

which are nicely satisfied by our coefficients up to machine double precision. This provides 



a non-trivial check of our numerical determination. Alternatively to Eq. 3.12 one could 
have used an expansion in spherical harmonics. The volume condition implies that the 
coefficient of Yoo(6,ip) should vanish. 
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We will now analyse the two-vortex dynamics that follows from our fifth order metric. 
This is governed by the Hamiltonian 



1 ■>. . vl 



H=T7^[P 2 e + zth) (3-14) 



The equations of motion are 



f(0,<p) \ sin 2 
f 



(3.15a) 



<P = (3.15b) 
j sin 9 



1 / ^ d f 9 cos 9 
f \ o9 v svn 6 



p e = -,[E^ + 2 P l^- (3.15c) 



Ptp = j— (3.15d) 

Where E is the energy of the system, that is conserved. These equations can be easily 
integrated numerically to obtain the trajectories in vortex moduli space. 

Notice, that up to third order, the metric is independent of (p. Hence, p v is a conserved 
quantity, and the problem is integrable. In this case, the equation for the trajectory in the 
moduli space is given by 



I s d9 

J0o c ;„2 njm 1 



00 Sill u \i -T2 ^r- a 

1 b A sin 2 



where b = ^fi=. 

An interesting mathematical question is whether the two vortex dynamics on the torus 
is integrable for all values of e. As we have seen, this is the case up to third order in the 
expansion. In the next subsection we will mention that it is also the case in the e — > 1 limit, 
since then we recover rotational invariance and angular momentum becomes a conserved 
quantity. Our metric calculations could be used to give a numerical and/or analytical 
answer to this question. We will not attempt to study that point here. Instead we will 
address the much simpler problem of obtaining and describing the Poincare maps obtained 
from the third order and fifth order metric. Results are displayed in Fig. [l] for two values 
of e. We see signs of how some invariant-tori seem to be destroyed, pointing to the non- 
integrability of the metric truncated to fifth order. This is not surprising and does not 
conflict with the conjectured integrability for all values of e. Our results, nonetheless, can 
be considered as a first step, which might help in pointing to the relevant region of phase 
space where one should look. 

3.2 Dynamics of vortices in the plane (e — > 1 limit) 

The dynamics of the Abelian Higgs model in M? is a well studied problem, both ana- 
lytically |2(^, ^IJ and numerically [22, pif| . Following the notation of Ref. |^TJ we might 



write 

„2 Q 77*2 



ds 2 = 8F 2 (2|n|) dudu (3.17) 
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(a) Poincare map for e — 0.4 and E — 60.0 for 
the first 3 orders. 



(b) Poincare map for e = 0.4 and E = 60.0 for 
the first 5 orders. 




E=200.0;e=0.4 



0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 

cos(O) 



(c) Detail of the Poincare map for e = 0.4 and (d) Detail of the Poincare map for e = 0.4 and 
E = 200.0 for the first 3 orders. E = 200.0 for the first 5 orders. 

Figure 1: Examples of Poincare maps for the dynamics near the Bradlow limit. 



Once the function F^u]) is obtained, the dynamics in the geodesic approximation is 
completely determined. The form of the function .F(2|u|) has been deduced analytically 
for the special case of asymptotically separated vortices [27], and has also been studied 



numerically several times [20, 21]. In this paper we will use our expansion method to 
obtain a precise determination of -F(2|u|) for all values of |u|. 

We will naturally assume that, if the torus is large enough compared both with the 
typical size of vortices, and with the distance between them, the dynamics will be accurately 
described by the one in M 2 . This means that, at least formally, we can calculate the desired 
function _F(2|u|) by taking the e — ► 1 limit while keeping u fixed. At the beginning of 
this section we showed how to express the metric in the relative torus position coordinate 
w. The formula (Eq. |3.8| ) involves f w times a known function arising from the change of 
variables. We will first study the behaviour of the latter factor in the appropriate limit 



(e — > 1 , u fixed). Restricting to the l\ = I2 = I case and recalling that w = Uy/2ir(l — e) 
we get: 

^jW a =4 - e)W (3.18) 

W 1 1 ;I (|^HI 2 + |^MI 2 ) 2 (i + ^o 2 ) 2 

where 3 Z 2 = -% M £$ 0) , and Z = -tf 3 (0)A? 2 (0). 

In summary, the function F 2 (2\u\) can be recovered as 



2vr 2 Z 2 



F 2 (2|u|) = lim (1 - e) 2 / tu (n v / 2vr(l - e),u v / 2ir(l - e), e) „ ^ 2 > | 2 (3.19) 



In practice, to obtain an approximation to -F 2 (2|w|) from our expansion we proceeded as 
follows. We computed the coefficients in the expansion in e of the metric up to 40th order, 
for a regular rectangular array of points in the complex plane of the variable w. The range 
of values is that of a complex torus. Due to symmetries (which we checked) it is enough 
to take values in one quadrant. We considered a grid of 40 x 40 points. Getting to order 
40 is feasible by implementing the iterative solution of the equations on a computer. This 
is possible using a Fourier expansion of the functions h and H and truncating to a finite 
number of Fourier coefficients. One can tune this finite number until the precision attained 
in the function / is of 14-16 decimal places. Due to the fast convergence of the Fourier 
coefficients, this procedure is not very costly. The computation in a conventional PC takes 
around 50 hours. 

Once / is determined at the grid points, in order to extract the values of i ?2 (2|u|) we 
proceeded as follows. For a given value of \u\ we take many different values of w from the 
grid and for each we compute the corresponding value of e. For that value we sum the 
series up to 40th order and multiply by the function associated to the change of variables 
(Eq. |3.8j ). Finally, for each value of \u\ we obtain we obtain a collection of numbers for 
different values of e and for different phases of u. A typical case is displayed in Fig. |2|, 
where we considered points having \u\ = 1 and situated along the diagonal Re(u) = Im(u). 
The value of F 2 (2| u\) is the limit for e approaching 1. Unfortunately, the larger values of e 
are more seriously affected by the truncation of the series to 40 orders. This is exemplified 
by the errors attached to the points of the figure. These are obtained by subtracting the 
result obtained for 35 orders to that of 40 orders. To obtain a more precise determination 
of the value of -F(2|u|) we made two corrections. First of all, we partially corrected for 
the truncation errors by estimating the contribution of the terms from the 41st on. The 
precise way in which this is done will be explained later. This additional contribution 
does not affect significantly points for which the error obtained by comparing the 35th and 
40th order is tiny. On the opposite extreme the contribution becomes unreliable when this 
error becomes very large. This occurs typically for values of e exceeding 0.95 (areas which 



3 For the rectangular symmetric case li = h, the numerical values of Zq and Z2 are 

Z = -1 - V2 

Z 2 = -1.18 91 73 93 79 71 08 15 
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are more than 20 times the critical area). We will therefore omit those points from our 
analysis. 

To extract a good determination of F 2 (x) from our data we face the problem of ex- 
trapolating to e = 1. This implies relating the results for a large torus with that of the 
plane. The torus case can be considered a solution of the plane with an infinite number of 
vortices - two per I x I cell. As the torus gets larger the replica vortices are increasingly 
far away, so that it is reasonable to use the analytic result for the case of a large number of 
asymptotically separated vortices to describe the approach. Using the result of Ref. j27]] we 
expect an additional contribution to the metric proportional to ^2 r Ko(D r ), where Kq is 
the modified Bessel function of the second kind and D r is the distance to the replica vortex 
r. In the limit of large torus sizes and fixed |u|, this predicts a contribution proportional 
to Kq(1). On the basis of the asymptotic behaviour of Bessel functions we decided to fit 
our data for each value of u to a formula of the form 



A-B(l- e)4e 



V87r/(l-e) 



(3.20) 



For our \u\ = 1 example, the result is displayed by the solid line in Fig|| and the fitted 
parameters are B = 145.0(5) and A = F 2 (2) = 0.8791(4). The error reflects systematics 
due to changing the range of the fit as well as the relative weighting of the points. For 
other phases of u we get compatible results. For example, for purely real or imaginary 
values of u we obtain F 2 (2) = 0.8800(6), consistent with rotational invariance. 



0.89 




0.95 



Figure 2: The points are the values of the left hand side of Eq. |3.19| for |u| = 1.0 as a function 
of e, computed from the first 40 orders of our expansion. Error bars represent the size of the last 5 



terms of the expansion. The solid line is a fit with Eq. 3.20 



Proceeding in the same way for several values of |it| we arrive at the points displayed 
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in Fig. ||. The errors are smaller than the size of the points in the figure. In all the range 
up to \u\ > 2 the error is smaller than 6 10 -4 and for \u\ > 0.5 the relative error is smaller 
than 1 part in 10 3 . 




Figure 3: The metric computed using 40 orders of the Bradlow parameter expansion (red points), 



versus the |it| — > approximation Eq. 3.25 (green dashed line) and the asymptotic form of the 



metric Eq. |3.2q (black solid line) . 



It is quite interesting to analyse the implications of the scaling limit for the coefficients 
of the e expansion. For that purpose it is convenient to write e = e~ s and replace the limit 
by 5 — ► 0. The condition that at large volumes the metric tends to that of the plane 
implies: 

F 2 (2\u\) = W %^ 2 f| 9 > | 2 f; f( N \uV2ri,uV2ri)e- NS (3.21) 

N=l 



6^0 (I + Z 2 



It is clear that as 5 tends to zero, more terms in the series become relevant. It is tempting 
to assume that the infinite sum tends to an integral over the variable y = N5. This would 
imply the following behaviour of the coefficients: 



(3.22) 



In this case the limiting function h(\u\^/y) would be related to F 2 (2\u\) as follows: 



F 2 (2|n| 



(i + zlY 

which is essentially a Laplace transform. 



~Vh{\u\Jy)dy 



(3.23) 
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With our results for up to N = 40 and the set of values of w spanning the 



aforementioned 40 x 40 grid, we have analysed the validity of Eq. 3.22. For that purpose 



we plot in Fig. || the coefficients satisfying 35 < N < 40 and all values of u in the range. 
The smoothness and thinness of the resulting curve, approximating h(x), is a test of our 
scaling hypothesis. Errors, signalled by the spread of the values, grow with x. For example, 
including all coefficients N > 20 would provide little changes in the thickness of the curve 
up to x ~ 5. Beyond, the thickness becomes a sizable fraction of the value. 

The function h(x) for small x, as determined from our results, fits nicely to a behaviour 
of the form 

h(x) = a-bx 4 + 0{x 5 ) (3.24) 

where a = 0.335(1) and b = 0.16(1). This expansion implies that for small values of the 
argument the metric function i ?2 (2|n|) behaves as: 



167T 2 Zf 



F 2 (2\u\ 



(a\u\ 2 - 2b\u\ 6 ) + 0{\u\ 7 ) 



(3.25) 



(1 + z%y 

The function on the left-hand side of the previous formula is displayed in Fig. ||[ It matches 
nicely to the behaviour of the previously determined points (in red). For large values of \u\ 
the data also matches with the prediction coming from the asymptotic behaviour at large 
distances (Ref. p7|| ) given by 



F 2 {2\u\) lu ^° 1 - 8V2K (A\u\) 



(3.26) 



where Kq{x) is the modified Bessel function of the second kind. 




Figure 4: We plot the left-hand side of Eq. 3.22 for N in the range 35-40, and fixed x — \u\^fy 
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To conclude we point out that our knowledge of h(x) enables one to estimate the error 
introduced by the truncation in the number of terms in the e expansion. For that one has 
simply to use the expression given in Eq. 3.23) with the integral restricted to y > —Nlog(e), 



where N is the maximum order computed. We have tested this estimation with our data 
points and iV = 40. As mentioned previously, it improves the results for not too large 
values of e. 

3.2.1 Scattering of vortices on the plane 

Once the function _F(2|u|) has been obtained, it is trivial to compute the dynamics of the 
vortices. The Lagrangian, when written in polar coordinates 2u = re ld , is given by: 

T = irF 2 (r)\f 2 + r 2 6 2 ] (3.27) 

Since the energy is conserved, and we have rotational invariance (pg is also conserved), the 
system is integrable, and we have 

r° dr . 

6 f = ir- (3.28) 

Jr f r v /^F2(2r)-l 

2 

where b 2 = is the impact parameter. In particular, if we call r m the minimum distance 
between one vortex and the centre of mass during the trajectory, the scattering angle is 
given by 

f°° dr 

e sc = 7r-2 (3.29) 
Jr m r J£ F 2( r }_i 

Using our numerical data of Fig. || we can calculate these integrals. For values of r > 5 
we will use the asymptotic form of the metric, and for r < 0.5 we can use our approxi- 



mation Eq. ( 3.25 ). For the intermediate values we will use the cubic spline interpolating 
polynomial. In the Fig. [5] we have some examples of trajectories for different values of the 
impact parameter. 

We can also plot the scattering angle versus the impact parameter. The results are 
shown in Fig. ||. Errors were estimated in the following way. We generate a random set of 
values of -F(2|it|) assuming a Gaussian distribution around the central value and standard 
deviation fixed by the error. For each realization we compute the cubic spline interpolation 
polynomial and using it and the small and large approximation functions we compute the 
scattering angle for each impact parameter. The final errors are estimated on the basis of 
1000 realizations. They are of the order of 1 part in 10 4 or smaller. 

4. Conclusions 

In this paper we have shown how the Bradlow parameter expansion can be used to compute 
the metric of the moduli space of solutions of the Abelian Higgs model on the Torus. This 



extends our previous work [16| where this expansion was introduced as a tool to construct 



the classical vortex solutions on the Torus. More precisely, the expansion is performed in 
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Figure 5: Scattering trajectories of two vortices. 
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Figure 6: Scattering Angle vs. impact parameter. 



the variable e = 1 — 2/, where / is the average magnetic field. To any finite order in the 
expansion the metric has a known form depending on a small number of real parameters. 
We give their analytic expression up to second order for arbitrary number of vortices q. The 



leading term, giving the Bradlow limit result, is proportional to the Fubini-Study metric. 
This coincides with results obtained for the two-sphere |3(|. The second order result 
depends on q 2 constants, which are given in the form of rapidly decreasing double infinite 
sums. Although increasingly complicated it is not too hard to go beyond this order. Rather 
than attempting this for the general vortex case, we have concentrated on the two-vortex 
case. We have given the form of the metric up to 5th order and the numerical value of their 
15 real coefficients with 14-16 significant digits (machine double precision). Furthermore, it 
is possible to go much beyond this order in the computation of the metric at specific points 
of the moduli space by means of a computer. Essentially, the implementation is based on 
the truncation of the infinite sums to finite ones. In this way we have gone up to 40th order 
with fairly limited computational effort. Furthermore, the fast convergence of the sums 
stills maintains the double machine precision accuracy in the determination of the metric. 
Having so many terms in the expansion has allowed us to extrapolate our results to the 
infinite area case (vortices on the plane) with quite impressive results. The resulting metric 
is consistent with all the previously known results, both analytical in certain limits and 
numerical. The result have fairly small and controlled errors, arising from the extrapolation 
to infinite volume and the truncation in the expansion. 

The computation of the metric allows the study of vortex scattering, known to be a 
good approximation up to relatively high vortex velocities. The presence of the torus can 
be used to mimic the behaviour of vortices in a dense environment. A natural question 
that poses itself is the possibility that two- vortex dynamics on the torus is integrable. This 
is actually the case for both small (close to critical) torus sizes as for infinitely large ones 
(the plane). We have analysed Poincare maps as a first step in this study. Furthermore, we 
emphasise that, as in our previous work, our method has been found to be a competitive 
tool to study the properties of vortices on the plane. It has important advantages over 
other methods: its semi-analytical nature and the possibility of dealing with an arbitrary 
number of vortices. In this work, for example, we have given the curves (with errors) 
that describe the functional dependence of the two- vortex scattering angle with the impact 
parameter. 

To conclude, we want to mention that the Bradlow expansion method is generalisable 
to many other cases and theories. For example, it can be applied to any other two- 
dimensional compact manifold. Furthermore, a similar expansion seems to apply in other 
cases where a Bradlow parameter can be introduced []l7|], such as abelian and non-abelian 
gauge-Higgs theories in Kahler manifolds of arbitrary dimensionality. The methodology 
can be extended beyond. For example, a related expansion, is known to exist in the study 
of self-dual configurations on the 4 dimensional torus [18]. Recently [35, 3C], it has also 



been applied to the study of abelian vortices in a non commutative torus. In all these 
theories, apart from the determination of the shape of the solutions and the form of the 
metric of the moduli space, other problems can be attacked. In particular, one can study 
the eigenfunctions and eigenvalues of the Dirac operator in a vortex background and also 
compute quantum corrections to the vortex. These problems are currently under study by 
the present authors. 
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A. Apendix: 

To lowest order we have that /i (0) = m(||c||) and # (0) = are constant. The symbol 

||c|| 2 stands for ^ |cj| 2 . Hence, the lowest order metric is given by 

The matrix inside parenthesis, which we will call V, is the orthogonal projector to the 
vector Cj. It gives precisely the Fubini-Study metric written in homogeneous coordinates. 
To the next order the metric, written in matrix notation, adopts the form: 

7r||c|p \ \\c\\ z I 

where 

Sij = 2 / ^XiXj^iXiXk)' ( A -3) 

and 

= J^I / d'xXiXkA-'ixiXj)' (A.4) 

The prime denotes the removal of the constant Fourier mode, making the inverse laplacian 
A -1 well-defined. 

The computation of the integrals can be explicitly carried over by using Fourier trans- 
forms. For that purpose we introduce the quantities: 



1 



,2 _ -2mni--A 



L*(ft) = - I d'xx&^Xi (A-5) 



T 2 



which are just the Fourier modes of the product XiXj- Now we remove the constant mode 
(n = 0) and denote the result as L' l3 {n). Then the integral entering in the definition of S 
and S is given by 

s m - ~\ / T2 ^x,A- W)' = £ ^^Q^ ( A - 6 ) 



where 

t(n) = — n 2 + -| 



(A.7) 
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and r = k/h- The explicit value of the coefficients depends on the choice of basis. For the 
choice given earlier Xs+i = ^Os defined in equation (A24) of our previous paper ||l6fl , the 
coefficients are given (as a particular case) of Eq. (A. 30) of that paper. The matrix L(n) 
becomes: 

tin) 

(A.8) 



((H) 

L(n) = r(n)e 2 



where the matrices T(n) are defined as 

T(n) = P ni Q n ''e 



. "1"2 
1 



(A.9) 



where P and Q are 't Hooft q x q matrices (Pij = <5ji+i, Q = diag(exp 2nis/q)). The 
matrices T(ft) are unitary and satisfy T^(n) = T(—n). From here one realises that, the 
coefficients SV,/y can be expressed in terms of q 2 constants K(s) as follows: 



ijkl 



E 



(A.10) 



where s is a 2 component vector with components Sj G {0, . . . ,q— 1}. From Eq. A. 6 
derives: 



one 



fcez 2 



e -?(s+gfc) 



da 



si/q 




(0\iqar) •& 



S2/q 




(0;iga/r)-l (A.ll) 



The terms in the sum decay exponentially with k 2 , so that a few terms in the infinite sums 
over k suffice to obtain a very high accuracy. 

The matrices T(s) are a basis of the set of q x q matrices. Hence, it is possible to 
express the metric tensor to this order as a linear combination of these matrices. The 
coefficients are functions of the following quadratic forms in Cj/||c||: 



P(s) 



cT(-s)c 



(A.12) 



These are not independent, and satisfy certain relations. In particular one has P (0) = 1 
and 



Y^\P(s)\ 2 = q - 1 

As an example, we give the explicit expression for the trace of g^: 
||c|| 2 Tr( 5 ( 2 )) = 4 £ JT(3)((, + 1)|P( S -)| 2 - 1) 



(A.13) 



(A.14) 



From here it is trivial to reproduce the q = 2 case given in section 3.1. Notice that 
(P(l, 0), P(0, 1), -P(l, 1)) is a 3-component unit vector whose expression in terms of spher- 
ical coordinates coincides with our definition of 9 and ip in Eq. For the symmetric 
h = h case we obtain 



/ 



(2) 



A 



IT 



(^(1,0)-K(1,1))(1-3|P(1,1) 



(A.15) 
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where P(l, 1) = cos 9. This reproduces the result given in Table |]. Notice that the result is 
proportional to the second Legendre polynomial in cos 9 as follows from the area condition. 

Going over to higher orders is tedious but presents essentially no new difficulties. We 
have obtained the analytic result to third order for the q = 2, r = 1 case, matching the 
value given in Table [|. 
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